Functionalization of Single Layer MoS 2 Honeycomb 

Structures 



C. Ataca 1 ^ and S. Ciraci*'* 

Department of Physics, Bilkent University, Ankara 06800, Turkey, and UNAM-Institute of 
Materials Science and Nanotechnology, Bilkent University, Ankara 06800, Turkey 

E-mail: ciraci@fen.bilkent.edu.tr 



*To whom correspondence should be addressed 

^Department of Physics, Bilkent University, Ankara 06800, Turkey 

^UNAM-Institute of Materials Science and Nanotechnology, Bilkent University, Ankara 06800, Turkey 
^UNAM-Institute of Materials Science and Nanotechnology, Bilkent University, Ankara 06800, Turkey 



1 



Abstract 

Based on the first-principles plane wave calculations, we studied the functionalization of 
the two-dimensional single layer M0S2 structure via adatom adsorption and vacancy defect 
creation. Minimum energy adsorption sites are determined for sixteen different adatoms, each 
gives rise to diverse properties. Bare, single layer M0S2, which is normally a nonmagnetic, 
direct band gap semiconductor, attains a net magnetic moment upon adsorption of specific 
transition metal atoms, as well as silicon and germanium atoms. The localized donor and ac- 
ceptor states in the band gap expand the utilization of M0S2 in nanoelectronics and spintronics. 
Specific adatoms, like C and O, attain significant excess charge upon adsorption to single layer 
M0S2 which may be useful for its tribological applications. Each MoS2-triple vacancy created 
in a single layer M0S2 gives rise to a net magnetic moment, while other vacancy defects related 
with Mo and S atoms do not influence the nonmagnetic ground state. Present results are also 
relevant for the surface of graphitic M0S2. 

Introduction 

The synthesis of molybdenum disulfide, M0S2, flakes have been subject of interest due to the 
unique symmetry of honeycomb structure and esoteric binding mechanism. 1-3 Recently, 2D sus- 
pended single layer M0S2 sheets (specified as IH-M0S2) have been produced.-^ Single layer 
M0S2 nanocrystals of ~ 30 A width were also synthesized on Au(lll) surface. 2 Liquid exfo- 
liation of nanosheets of M0S2 and other transition metal dichalcogenides, such as MoSe2, WS2, 
MoTe2, NbSe2, NiTe2, are reported. 6 IH-M0S2 is made of hexagons with Mo and two S atoms 
situated at alternating corners. Recently, this material with its nanoribbons has been an active field 
of theoretical studies.-"— 

M0S2 has been used in ultra low friction studies. Formerly, the surface of 3D graphitic M0S2 is 
used for these experiments to achieve a friction coefficient of 10 -3 along the sliding direction in its 
basal planes. 13 Recently, experiments-^^ using friction force microscopy and theoretical calcula- 
tions — on atomically thin sheets of M0S2 are confirmed the superlubricity of hexagonal M0S2 
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structure. In addition to superlubricity, other areas, where M0S2 appears to be a potential material 
for various technological applications, are hydrogen production,-^*^ hydrodesulfurization catalyst 
used for removing sulfur compounds from oil, 20-26 solar cells, 27 photocatalysis. 28 Photolumines- 
cence is measured in monolayer M0S2 which is absent in graphitic structure. 29 Most recently, a 
transistor fabricated from single layer M0S2 has indicated the superior features of this 2D material 
over graphene. 30 Studies to date suggests that M0S2 can be promising for optoelectronic devices, 
solar cells, and LEDs. 

Various properties revealed earlier for 3D graphitic M0S2 and nowadays for 2D single layer 1H- 
M0S2 have made the functionalization of these structures by adatoms or vacancy defects a current 
issue. Ab-initio calculations were carried for adsorption of thiophene 31 on catalytically active 
surface of M0S2 crystal, which is used in hydrodesulfurization process. Similarly, the adsorption of 
aromatic (thiophene, benzothiophene, benzene, naphthalene, pyridine, quinoline) and conjugated 
compounds (butadiene) on the basal plane (0001) of M0S2 were studied. 32 Huang and Cho 33 
investigated the adsorption of CO on M0S2 surface and calculated the relative energies of different 
reaction paths. Implementing local magnetism through defects or impurities has been the focus of 
first-principles calculations. Fuhr et al. 34 found that S vacancy defect on M0S2 surface, as well 
as substitutional doping of Pd and Au do not induce magnetic moment, whereas Fe and V induce 
magnetic moments when substituted with S atoms at M0S2 surface. Adsorption and substitutional 
doping of Nb atoms on M0S2 surface were also investigated. 35 Since the magnetism based on 
s/?-orbitals yields long-range exchange coupling interactions as compared to d— and /— orbitals, 
implementing magnetic properties to M0S2 monolayer through nonmagnetic adatoms have been 
also considered. Theoretical studies concluded that H, B, C, N and O can significantly modify the 
magnetic and electronic properties of this material. 36 

In this paper we present a systematic study to reveal the general effects of the adsorption of 
selected adatoms and vacancy defect on the electronic and magnetic properties of perfect 1H- 
M0S2. Three issues are of particular interest: (i) How can the surface charge of the single layer of 
M0S2 be enhanced through adatoms? (ii) How can the spin polarization and magnetic moment be 
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induced in nonmagnetic IH-M0S2? (iii) How can the electronic and magnetic states be modified 
through adatoms and vacancy induced localized gap states? We investigated 16 individual adatoms 
(such as C, Co, Cr, Fe, Ge, Mn, Mo, Ni, O, Pt, S, Sc, Si, Ti, V and W) and 5 different types 
of vacancy defects (such as Mo-, S-vacancy, S2-, MoS-divacancy and M0S2 triple vacancy) in 
view the above three issues and investigated their general features. Owing to the weak interlayer 
interaction present predictions concerning functionalization of IH-M0S2 honeycomb structure are 
expected to be relevant also for the M0S2 sheets comprising a few layers, as well as the surface of 
graphitic M0S2. 

Method 

Our results are based on first-principles plane wave calculations within density functional theory 
(DFT) using projector augmented wave (PAW) potentials. 37 The exchange correlation potential 
is approximated by generalized gradient approximation (GGA) using PW91 38 functional both for 
spin-polarized and spin-unpolarized cases. While all discussions in the paper are based on the re- 
sults obtained within GGA using PAW potential, calculations within local density approximation 39 
(LDA) using PAW are also performed for specific cases for the purpose of comparison. All struc- 
tures are treated within the periodic boundary conditions using supercell geometry. Kinetic energy 
cutoff, Brillouin zone (BZ) sampling are determined after extensive convergence analysis. A large 
spacing of ^ 10 A between 2D single layers of M0S2 is taken to hinder the coupling between them. 
A plane- wave basis set with kinetic energy cutoff of 600 eV is used. In the self-consistent field 
potential and total energy calculations, BZ is sampled by special k-points.— The numbers of these 
k-points are (35x35x1) for the unitcell and (5x5x1) for the adatom adsorption in (4x4) supercell of 
IH-M0S2. All atomic positions and lattice constants are optimized by using the conjugate gradient 
method, where the total energy and atomic forces are minimized. The convergence for energy is 
chosen as 10 -5 eV between two consecutive steps, and the maximum Hellmann-Feynman forces 
acting on each atom is less than 0.02 eV/A upon ionic relaxation. The pressure in the unit cell is 
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kept below 1 kBar. Numerical calculations have been performed by using VASP. — & 

Since DFT within GGA underestimates the band gap, frequency-dependent GWo calcula- 
tions 43 are carried out to correct the band gaps. Screened Coulomb potential, W, is kept fixed 
to initial DFT value Wq and Green's function, G, is iterated four times. Various tests regarding 
vacuum separation, kinetic energy cut-off energy, number of bands, k-points and grid points are 
made. Final results of GWo corrections are obtained using (12x12x1) k-points in BZ, 400 eV 
cut-off potential, 192 bands and 64 grid points for IH-M0S2. We were not able to apply GWo 
corrections for adatoms and vacancy due to large number of atoms in the supercells. 

(a) 

Top View Side View 
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Figure 1: (a) Top and side views of atomic structure of 2D IH-M0S2 with hexagonal lattice. The 
hexagonal unit cell with lattice constants \a\ = \b\ is delineated by thin solid lines. Honeycomb 
structure consisting of Mo and S2 atoms located at the corners is highlighted by dotted hexagons, 
(b) Contour plots of charge density, p (see text for definition) in a vertical plane passing through 
Mo-S bonds. Arrows indicate the increasing value of charge density, (c) Isosurface plot of differ- 
ence charge density, Ap (see text for definition). Isosurface value is taken as 0.006 electrons/A 3 . 

(d) Energy band structure of IH-M0S2 calculated by GGA+PAW using optimized structure. Zero 
of energy is set to the Fermi level indicated by dash-dotted line. The gap between valence and con- 
duction band is shaded; GWo corrected valence and conduction bands are shown by filled circles. 

(e) Total density of states (TDOS) and orbital projected density of states (PDOS) for Mo and S. 
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Table 1: Calculated values of single layer IH-M0S2: Lattice constants, \a\ = \b\; bond lengths, 
dMo-s and dss\ S-Mo-S bond angle, 0(S-Mo-S); cohesive energy per M0S2 unit, Ec; direct 
energy band gap, E g ; band gap with GWq correction, E g+ cw ; photoelectric threshold, <E>. 



Method 


a (A) 


dMo-S( A ) 




fl(S-Mo-S) 


E c (eV) 


E g (eV) 


E g+GW {eV) 


4>(eV) 


GGA+PAW 


3.20 


2.42 


3.13 


80.69° 


15.55 


1.58 


2.50 


5.96 


LDA+PAW 


3.11 


2.37 


3.11 


81.62° 


19.05 


1.87 


2.57 




Experiment 


3.27, 44 3.20 45 
















Theory 


3.13 7 


2.41 7 




82.31° 7 


15.6, 7 14.99 9 


1.69, 7 1.8, 9 1.78, 8 1.64 10 







Two dimensional honeycomb structure of MoS 2 

To reveal the effects of adatoms and vacancy defects, we first present our results related with 
atomic, electronic and magnetic properties of IH-M0S2, which are obtained using the same calcu- 
lation parameters. In |[figure][l][]l| (a) we schematically illustrate the hexagonal crystal lattice, re- 
sulting in a honeycomb structure. The contour plots of total charge density, p(r) = Yf cc |*F(k,r) | 2 , 
and difference charge density isosurfaces, Ap(r) (which is obtained by subtracting the free atom 
charge densities of Mo and S atoms from the total charge density of IH-M0S2) are also shown 
[figure] [!][]! | (b) and (c), respectively. There are ambiguities in determining the charge trans- 
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fer; calculated excess charges may depend on the method used. While different methods result in 
different values for charge transfer (or excess charge on the atoms), its direction can be predicted 
unambiguously. Based on Mulliken— ^ analysis the excess charge on each S atom and depletion 
of electrons on each Mo atom is calculated to be 0.205 electrons and 0.410 electrons, respectively. 
Bader analysis 48 yields depletion of 1.2 electrons on Mo atom and 0.6 excess electrons on each 
S atom. Electronic charge transferred from Mo to S atoms gives rise to an excess charge around 
S atoms. Consequently, IH-M0S2 can be viewed as a positively charged Mo planes sandwiched 
between two negatively charged planes of S atoms as presented in |[figure][l][]l| (c). The repulsive 
Coulomb force between negatively charged S atoms facing each other in adjacent M0S2 layers 
weakens the interlayer interaction in two individual IH-M0S2 (or bilayer of M0S2), as well as 
in graphitic M0S2. In fact, the interaction between two adjacent M0S2 layers of graphitic M0S2 
is repulsive even for large separation except a slight attractive range (of only 7 meV without van 
der Waals interaction, but ^170 meV including van der Waals interaction) around the equilibrium 
separation. This is the main reason why M0S2 flakes can be used to lower the friction coefficient in 



sliding friction. The cohesive energy of the optimized IH-M0S2 structure is found to be 15.55 eV 
per cell. For the known results, LDA+PAW calculations yield stronger binding and hence larger 
cohesive energy (£c=19.05 eV). 

Electronic structure of 2D suspended single layer IH-M0S2 and corresponding TDOS and 



PDOS calculated within GGA+PAW are presented in | [figure] [l][]l| (d) and (e). The states at the 
band edges, as well as at lower lying conduction and higher lying valence bands originate from the 
hybridization of Mo-4d and S-3p orbitals. The direct band gap is £^=1.58 eV. The band gap cal- 
culated within LDA+PAW is E g =1.87 eV, which is in good agreement with the band gap measured 
using complementary techniques of optical absorption, photoluminescence and photoconductivity 
as 1.90 eV. 5 Earlier Bollinger et al. 10 predicted the band gap to be 1.64 eV using GGA+ultrasoft 
pseudopotential. However, Li and Galli 12 with similar computational method calculated the band 
gap as 1.80 eV. Li et al. 1 calculated the direct band gap 1.69 eV within GGA+PAW using rela- 
tively smaller energy cut-off. Mendez et al. 9 found the direct band gap within LDA as E g =l.% 
eV using local basis set. 49 Lebegue and Eriksson 8 carried out LDA+PAW calculations using ex- 
perimental lattice constants and found the band gap to be 1.78 eV. The band gap calculated in 
the present work is in fair agreement with previous studies.-^ However, we show that band 
gaps determined in earlier studies increase ~1 eV upon GWo correction. The band gap calculated 
within GGA(LDA)+PAW is corrected using self-energy method GW to be 2.50(2.57) eV. The 
corrected band gap is ^0.6-0.7 eV larger than the value measured experimentally. 5 The situation 
with graphitic M0S2, which consists of the stacking of IH-M0S2 layers is, however, different. 
The indirect band gap of graphitic MoS 2 calculated within GGA(LDA)+PAW is 0.85(0.72) eV 
and is corrected to be 1.44 (1.28) eV. In particular LDA+PAW band gap corrected by GWo is in 
good agreement with the experimental value. 50 Since LDA/GGA is designed to describe systems 
with slowly varying electron density and may fail to model localized d-orbitals, 51 more accurate 
band gap calculations can be carried out using generalized Kohn-Sham scheme, screened nonlocal 
exchange functional, HSE. 52 Starting from nonlocal charge density and wavef unctions, we calcu- 
lated HSE and Go Wo— corrected electronic band structure of M0S2 as a direct band gap of 2.23 



7 



and 2.78 eV, respectively. Surprisingly, the band gap of 2D IH-M0S2 is overestimated by GWo 
correction. Here we note that the band gap of 2D fluorographene CF is also overestimated upon 
GWo correction. 53 In stable] [1][]1| we list all the calculated structural parameters (including lattice 
constants, as well as internal parameters), cohesive energy, direct band gap together with GWo 
correction and photoelectric threshold of 2D IH-M0S2. 

Finally, we emphasize two important dimensionality effects related with the electronic structure 
of M0S2, which have important consequences like photoluminescence: (i) By going from graphitic 
M0S2 to single layer IH-M0S2 the energy band gap changes from indirect to direct, (ii) The 
minimum gap increases by ^0.6 eV. The transformation from indirect to direct gap is related with 
the orbital composition of states at the edges of conduction and valence bands.— ^ The widening 
of the band gap in 2D occurs due to the absence of S-p z orbital interaction between adjacent M0S2 
layers and appears as the manifestation of the quantum confinement in the direction perpendicular 
to the MoS 2 layer. 5 



Functionalization by adatom adsorption 
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Figure 2: Top and side views are the schematic representation of possible adsorption geometries of 
adatoms obtained after the structure optimization. Adatoms, host Mo and S atoms are represented 
by red (large-dark), purple (medium-gray) and yellow (small-light) balls, respectively. Side views 
clarify the heights of adatoms from Mo and S atomic planes. Different adsorption sites are specified 
below each entry as 'Mo(S)-(J', where Mo(S) indicates that adatoms are placed initially (before 
structure optimization) to Mo(S) plane. In Mo-1 and Mo-2 geometries the adatoms are in and 
slightly above the Mo-layer. S-l, S-2,.. S-4 positions are associated with the S-layer. The adatoms 
adsorbed at each site are given at the lower right hand side of each entry. 
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Adsorption of adatoms is widely used and an efficient way to provide new functionalities to 
structures in nanoscale applications.—"— Among 16 different adatoms, namely, C, Co, Cr, Fe, Ge, 
Mn, Mo, Ni, O, Pt, S, Sc, Si, Ti, V and W, we were interested in 10 transition metal elements to 
determine ones which can induce magnetization in nonmagnetic IH-M0S2. Group 4A elements, 
i.e. C, Si and Ge are considered, since these atoms can form stable, planar or buckled honeycomb 
structures.—"— Carbon was of particular interest for the fabrication of graphene+MoS2 complex or 
nanomeshes. Mo and S being host atoms in IH-M0S2, they can exist as residues. Whether M0S2 
can be oxidized is an important issue, which may limit future applications. The first question one 
has to address is, of course, whether these adatoms can form strong and stable bonding with 1H- 
M0S2. The equilibrium adsorption sites of these 16 atoms are determined by first placing them 
to one of four different adsorption sites and subsequently by optimizing whole structure. Four 
possible adsorption sites considered initially for each adatom before the structure optimization 
are: (i) Hollow site slightly above the center of hexagon at Mo atomic plane, (ii) Hollow site 
above the center of hexagon at S-plane. (iii) Top of Mo atom, (iv) Top of S atom. Here the 
bridge site above the Mo-S bond is very similar to adsorption of adatom at the top of Mo atom. 
In order to avoid adatom-adatom coupling, the (4x4) supercell of IH-M0S2 is used, where the 
nearest adatom-adatom distance is ~ 12.8 A. Thus, the results can be related to the adsorption of 
an isolated adatom. The minimum energy positions of adatoms after the optimization process are 



described in | [figure] [2] []2[ We note that C, Co, Fe, O, Ti, V and W have two different binding sites 
(one having higher binding energy, E^); each site leading to dramatically different electronic and 
magnetic structure. Six distinct adsorption sites described in | [figure] [2] []2[ namely Mo-1, Mo-2, 
S-l, S-2, S-3 and S-4 are distinguished after structure optimization. 

All relevant data obtained from our calculations of adatoms adsorbed to IH-M0S2 are presented 

in 



[table] [2] Q21 The height of the adatom from the Mo- or S -planes are calculated relative to the 



average heights of Mo- and S- atoms in the corresponding planes. The binding energy, is 
calculated as — E a ^ + Em s 2 ~ E a d+MoS 2 • Here, E a ^ is the ground state energy of free adatom 
calculated in the same supercell with the same parameters; Em s 2 * s the total energy of (4x4) 
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supercell of IH-M0S2 and E ac i+ Mo s 2 is the optimized total energy of adatom+(4x4) supercell of 



IH-M0S2. Among different adatoms studied here, Cr has the weakest binding energy =1.08 
eV); W has the strongest binding (Ejy =4.93 eV) and creates a local reconstruction on S-layer. 
The excess charge 48 of the adatom p*, is obtained by subtracting the calculated charge at the 
adatom, pA from the valence charge of the adatom Z^, namely p* — Za — Pa- Accordingly p* < 
implies excess electron at the adatom site. The magnetic moments are obtained by carrying out 
spin-polarized calculations without assigning initial atomic magnetic moments on adatoms. Upon 
relaxation the final magnetic moment is obtained. 

Table 2: Calculated values for the properties of 16 adatoms adsorbed on IH-M0S2. For 
specific adatoms, the first and second lines are associated with the adsorption to the Mo- 
layer and S-layer site, respectively. Other adatoms have only positive binding energy when 
adsorbed to the S-layer site. The adsorption sites of adatoms are described in |[figure][2][]2 
Hmo? the height of the adatom from Mo layer; hs, the height of the adatom from the nearest 
S-layer; dMo> the adatom-nearest Mo distance; ds, the adatom-nearest S distance; adatom 
binding energy; magnetic moment per supercell in Bohr magneton jis; p% excess charge 
on the adatom (where negative sign indicates excess electrons); photoelectric threshold 
(work function); P, dipole moment calculated in the direction normal to IH-M0S2 surface. 
E[, energies of localized states induced by adatoms. Localized states are measured from the 
top of the valence bands in eV. The occupied ones are indicated by bold numerals and their 
spin alignments are denoted by either t or |. States without the arrow sign indicating of spin 
alignment are nonmagnetic. 





Atom 


Site 


Hlo 


hs 


dMo 


d S 


E b 


Mr 


P* 


3> 


P 


Ej 






(A) 


(A) 


(A) 


(A) 


(eV) 


(Mb) 


(e) 


(eV) 


(ex A) 


t : Spin-up, | : Spin-down States 






Mo-1 


-0.01 


1.56 


2.04 


2.44 


3.28 


NM 


-0.65 


5.81 


0.00 


0.14, 0.19,2.18,4.85 


C 


S-3 


1.58 


-0.05 


2.07 


1.80 


2.69 


2.00 


-0.58 


5.74 


-0.10 


0.26 1, 0.30 1, 1.07 |, 1.13 t, 1-55 |, 1.60 | 


Co 


Mo-2 


0.61 


0.99 


2.22 


2.19 


0.96 


1.00 


0.15 


5.57 


-0.07 


0.61 1, 0.66 t, 0.82 1, 0.87 |, 0.89 1, 1.13 I 


S-3 


2.52 


0.94 


2.56 


2.11 


2.92 


1.00 


0.44 


5.17 


-0.38 


0.29 |, 0.51 1, 0.57 |, 0.62 |, 0.99 |, 1.00 | 


Cr 


S-3 


2.79 


1.21 


2.77 


2.26 


1.08 


4.00 


0.91 


4.25 


-0.75 


0.19 |, 0.22 |, 0.59 |, 0.64 f, 0.68 |, 1.60 |, 1.73 | 


Fe 


Mo-2 


0.31 


1.29 


2.21 


2.26 


0.39 


2.00 


0.34 


5.07 


-0.04 


0.32 |, 0.48 |, 0.87 |, 1.01 i, 1.15 |, 1.32 | 


S-3 


2.52 


0.93 


2.53 


2.14 


2.42 


2.00 


0.59 


4.99 


-0.46 


-0.03 t, 0.18 |, 0.21 1, 0.22 1, 0.91 1, 0.93 I, 1.3 H 


Ge 


S-4 


3.83 


2.26 


4.28 


2.30 


1.18 


2.00 


0.39 


4.78 


-0.27 


1.10 |, l.H t, 1-72 |, 1.75| 


Mn 


S-1 


2.51 


0.93 


3.06 


2.15 


1.37 


3.00 


0.81 


4.75 


-0.44 


0.16 |, 0.17 t, 0.49 f, 0.50 f, 1.30 1, 1.98 I, 2.05 4, 


Mo 


S-3 


2.89 


1.31 


2.84 


2.30 


1.43 


4.00 


0.81 


4.54 


-0.76 


0.48 |, 0.79 t, 0.87 |, 0.91 f, 1.16 1, 1.42 |, 1.54 I 


Ni 


S-3 


2.52 


0.94 


2.56 


2.12 


3.65 


NM 


0.36 


5.58 


-0.31 


0.12, 0.16, 0.21, 0.50, 0.51, 0.57 


O 


Mo-1 


-0.01 


1.57 


1.95 


2.56 


2.24 


NM 


-1.11 


5.72 


0.00 


-6.36, -6.34, -5.66, 0.32 


S-4 


2.94 


1.38 


3.51 


1.49 


3.99 


NM 


-0.91 


5.96 


0.35 


-5.63, -1.16, -0.90 


Pt 


S-3 


2.78 


1.21 


2.76 


2.31 


2.92 


NM 


0.08 


5.66 


-0.26 


0.26, 0.27, 0.37, 1.76 


S 


S-4 


3.48 


1.91 


3.96 


1.94 


2.17 


NM 


-0.11 


5.96 


0.23 


0.00, 0.06, 0.11 


Sc 


S-1 


2.84 


1.25 


3.30 


2.29 


2.63 


1.00 


1.45 


4.31 


-1.18 


1.28 |, 1.31 1, 1-35 I, 1.49 |, 2.25 I, 2.39 t 


Si 


S-4 


3.73 


2.16 


4.18 


2.17 


1.39 


2.00 


0.58 


4.94 


-0.14 


0.98 |, 0.99 1, 1-65 |, 1.68 4, 


Ti 


Mo-1 


0.00 


1.58 


2.31 


2.42 


1.23 


NM 


1.14 


5.71 


0.00 


0.26, 0.34, 0.36, 1.74, 1.77 


S-3 


2.95 


1.38 


2.99 


2.32 


2.62 


4.00 


1.16 


4.32 


-1.12 


1.18 |, 1.24 |, 1.35 |, 1.37 |, 1-54 f, 1-77 ^ 1.81 1 




Mo-1 


0.03 


1.55 


2.95 


2.34 


1.25 


1.00 


1.05 


4.61 


0.00 


0.12 |, 0.16 |, 0.17 t, 0.21 1, 1.54 1, 1.67 |, 1.69 4, 


V 


S-2 


1.87 


0.19 


2.56 


2.07 


2.76 


1.00 


1.05 


5.18 


0.04 


0.07 t, 0.18 |, 0.84 1, 1-23 |, 1.45 |, 1-54 | 


W 


Mo-1 


0.01 


1.58 


2.41 


2.42 


1.18 


2.00 


1.05 


4.66 


0.01 


0.16 |, 0.21 1, 1.41 1, 1.49 f, 1.61 ^, 1.81 4, 


S-2 


1.87 


0.16 


2.61 


2.15 


4.93 


NM 


0.85 


5.58 


0.12 


0.13, 0.19, 0.58, 1.57, 1.75 





Since the adatom-adatom interaction is hindered by a large separation between them, the 
adatoms presented in | [table] [2]G2| give rise to localized electronic states in the band gap and 
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resonant states in band continua and hence modify the electronic properties of IH-M0S2. In 



[figure] [3] [] 3, localized states of O, Ti, Cr and Ge together with band decomposed charge den- 



sity isosurfaces are presented. These atoms are specifically selected, since they are representatives 
of some of the adatoms presented | [table] [2] [] 2[ Oxygen adatom when placed on the S plane is ad- 
sorbed on top of S atom (S-4 site) with a binding energy of 3.99 eV. This site is in agreement with 
the results of He et al 36 However, we predict also a local minimum at Mo-1 site with relatively 
smaller binding energy (Et, =2.44 eV). Adsorbed O is nonmagnetic at both sites. Oxygen adatom 
having the highest electronegativity and highest negative excess charge among all other adatoms 
have localized states in the valence band. When O is adsorbed at Mo-1 site, the only localized state 
occurring in the band gap is filled and originate from the combination of 0-p z orbital with the 
/7-orbitals of nearest S atoms. Sulfur being in the same group with O displays similar electronic 
properties and have localized states in the band gap just above the valence band originated from its 
p x and p y orbitals. 



Titanium can be adsorbed to both Mo and S layers. In | [figure] [3] []3| localized states of both 
adsorption geometries are indicated. When Ti adsorbed to Mo-1 site, nearest d z i orbitals of Mo 
atoms play an important role in the binding mechanism. Charge density isosurface of the state 



numbered as '1' in [figure] [3] [] 3 is dominated from Ti's d z i, whereas isosurfaces '2' and '3' are 



dominated by d xy and d x 2_ y 2 of Ti which are parallel to Mo-layer. As for Ti adsorbed at S-3 site, 
Ti-d xy and d xz orbitals states dominate states numbered as ' V and '2' but the state numbered as '3' 
originates mainly from Ti-d z 2 orbital. Cr and Mo adatoms at S-3 site on IH-M0S2 surface have 
four d-states singly occupied, but d xy state is doubly occupied. Cr-s and Mo-s orbitals are vacant 
in both adatoms. Hence Cr and Mo has local magnetic moment of ji — 4/Ib- As for W adatom at 
S-2 site, its localized magnetic moment of /1 = 2/i# occurs from the spin-polarized d-orbitals of 
W together with the /?-orbitals of nearest S-atom. 

Carbon adatom is of particular interest. Previously, He et al 36 found that the lowest energy 
adsorption position of C adatom occurs in the plane of sulfur atoms on top of Mo atom (namely S-3 
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Figure 3: Schematic diagram of the relevant energy levels (or bands) of single adatom (O, Ti, 
Cr and Ge) adsorbed to each (4x4) supercell of IH-M0S2. The grey (light) shaded region in 
the background is the valence and conduction band continua. For nonmagnetic case, red (dark) 
bands are contributed more than 50% by adatoms orbitals, For magnetic case, spin-up and spin- 
down bands are shown by red (dark) and brown (light) lines, respectively. Solid bands indicate 
that the contribution of adatom to the band is more than 50%. In the lower part of each panel 
the adsorption site is indicated by the labeling of | [figure] [2] []2[ Charge density isosurfaces of 



adatom states specified by numerals are shown below. The isosurface value is taken as 2 x 10 5 
electrons/A 3 . 
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site shown in | [figure] [2] []2| ). Here, we found the same site with a binding energy of £^=2.69 eV. 
However, C adatom can go over a shallow barrier to a more stable adsorption position, namely Mo- 
1 site with £^=3.28 eV at the center of three Mo atoms below the plane of S atoms. In | [figure] [4] []4 



we present a detailed analysis of these two adsorption sites. For Mo-1 site, sp 2 hybridized orbitals 
of C adatom and d xy and d x 2_ y 2 orbitals of surrounding 3 Mo atoms form three bonds, but the 
system remains nonmagnetic. In the case of S-3 site, sp 2 hybridized orbitals of C adatom form 
three bonds with the orbitals of surrounding three S atoms of IH-M0S2. The interaction of p z 
orbital of C adatom with the underlying Mo d z 2 orbitals results in a significant charge transfer from 
Mo to C adatom. This creates a local magnetic moment of ji=2 jiB on Mo atom. The charge 
density contour plots in | [figure] [4] []4| (b) and (c) in a perpendicular and lateral planes show the 
bonding configuration. Recently a peculiar growth process, where carbon adatoms adsorbed to 
graphene readily diffuse above room temperature and nucleate segments of linear carbon chains 
(CACs) attached on graphene is predicted. 64 The diffusion barrier they calculate on graphene is 
reported as 0.37 eV. The energy difference of C adatom between different adsorption sites, namely 
S-3 and S-4, is 0.45 eV and comparable to that of C atom on graphene. Here we went beyond the 
adsorption of individual atoms and addressed the question whether carbon adatoms can also lead to 
the formation of CACs on IH-M0S2. Different from the single C adatom adsorption on IH-M0S2, 
CACs (starting from C2) prefer to adsorb on the top of S atoms as shown in | [figure] [4] []4| (d) and 
(e). When attached to IH-M0S2, CACs favor cumulene like structures with slightly alternating 
bonds and with a different kind of even-odd disparity. Therefore the top-site geometry occurs due 
to the double sp+ ^-bonding within cumulene. Similar to the case in graphene, the binding energy 
of CACs having even number of C atoms are greater than those having odd number of C atoms. 
Even though CAC on graphene are nonmagnetic structures, odd number CACs on IH-M0S2 have 
a magnetic moment of 2 /i# arising from C atoms at the edges of CAC. 

Since the atomic radii of other group 4A elements, Si and Ge, are larger than that of C adatom, 
their adsorption geometry differs from the adsorption site of C at the S-3 site. Si and Ge can 
only be adsorbed on the S -plane at the S-4 site and attain a local magnetic moment of /1 = 2 
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This is a crucial result, since magnetic properties to M0S2 monolayer are implemented through 
nonmagnetic adatoms. 



The excess charge on the adatom (which is denoted as p* in | [table] [2]"[]2| ), as well as the 
position of the highest localized state in the gap are correlated with its electronegativity. For 
example, among C, Si and Ge, C has the highest electronegativity. While C has negative excess 
charge, Si and Ge becomes positively charged when adsorbed. Calculated dipole moment on 
the system of C adatom adsorbed to S-layer is diminished due to the local reconstruction. While 
adsorbed C, O and S adatoms have excess electrons, all other adatoms in stable] [2] []2| are positively 
charged. Vanadium and W adsorbed at S-2 site have significant positive charge, but both induce 
minute electric dipole moment due to the local reconstruction. Finally, it should be noted that 
excess charging of surfaces of IH-M0S2 by a higher level of coverage of specific adatoms can 
improve its functionalities in tribology. 



Vacancy Defects 

It is known that the vacancy defect in 2D graphene, 65-69 graphene nanoribbons,—^ 2D graphane 72 
and graphane nanoribbons 73 give rise to crucial changes in the electronic and magnetic struc- 
ture. According to Lieb's theorem, 74 the net magnetic moment of the structure occurs due to 
the difference in the number of atoms belonging to different sublattices A and B, namely ji = 
(Na —Nb)lib- While magnetic moments calculated using DFT for vacancies in 2D graphene and 
armchair graphene nanoribbons generally confirmed Lieb's theorem, results are diversified for va- 
cancies in zigzag graphene nanoribbons-^^ due to the coupling with magnetic edge states. In the 
case of IH-M0S2, Lieb's theorem cannot be directly applicable, even if alternating Mo and S units 
form a honeycomb structure. We investigated five different types of vacancy defects, namely Mo-, 
S-vacancy, MoS-, S2-divacancy and MoS2-triple vacancy, which are formed in 2D IH-M0S2. All 
structures are optimized upon the creation of a particular type of vacancy. Vacancy energies, £y, 
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Figure 4: Analysis of bonding configuration of C adatoms on IH-M0S2. (a) Geometries of single 
C adatom adsorbed in the Mo-plane (left) and in the S-plane (right). Adatoms, host Mo and S 
atoms are represented by red (medium-dark), purple (large-gray) and yellow (small-light) balls, 
respectively, (b) Contour plots of total charge density of plane passing through atoms and bonds 
highlighted (not shaded) in (a), (c) Contour plots of total charge density on the horizontal plane 
passing through Mo-C and S-C bonds parallel to IH-M0S2. (d) Adsorption geometry and ener- 
getics of C2 and C3 on IH-M0S2. (e) Contour plots of total charge density on the vertical plane 
passing through atoms and bonds emphasized (not shaded) in (d). The direction of the arrows 
indicate the increasing charge density. 



are calculated by subtracting the total energy of the perfect structure (without vacancy) from the 
sum of the total energy of a structure having a particular type of vacancy and the total energy(ies) 
of missing atoms in the vacancy defect. Here all structures are optimized in their ground states 
(whether magnetic or nonmagnetic). Positive Ey indicates that the formation of vacancy defect 
is an endothermic process. In | [table] [3 ] [] 3 1 vacancy energies as defined above and their magnetic 
ground states are presented. We note that the equilibrium concentrations of vacancies are usually 
very low owing their high formation energies. Nonetheless, new techniques have been developed 
to create defects, which makes also the generation of nanomeshes vacancies possible.—"— 

Table 3: Calculated vacancy energies Ey (in eV), magnetic moments \i (in Hb) of five different 
types of vacancy defects, Mo, MoS, MoS 2 , S, S 2 in (7x7) supercell of 1H-MoS 2 . NM stands 
for nonmagnetic state with net /i=0 £/s denote the energies of localized states in the band 
gap measured from the top of the valence bands (in eV). The occupied ones are indicated 
by bold numerals and their spin alignments are denoted by either *f or |. States without the 
indication of spin alignment are nonmagnetic. 
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We carried out calculations on vacancy defects of Mo, S and S 2 , in the (2x2), (4x4) and (7x7) 
supercells of 1H-MoS 2 structure. Even if relatively larger sizes of supercells are not affected, the 
size of supercell is contracted upon creating a vacancy defect in the (2x2) supercell. In addition to 
these, we also consider MoS-divacancy and MoS 2 -triple vacancies in (7x7) supercells. For all these 
supercell sizes including the (7x7) supercell, we have carried out calculations with many different 
initial magnetic moment configurations, but all of the vacancy defects, except MoS 2 triple vacancy, 
the magnetic moment on any of the atoms in the supercell vanished. To justify these results, we also 
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repeated the calculations using LDA+PAW, which also results in a magnetic state for MoS2-triple 
vacancy having a net magnetic moment of \i? — 2/i#. We seek the origin of magnetic states in the 
charge transfer between Mo and S and hence examined charge densities around the vacancy defects 
using Bader analysis. 48 In case of S and S2 vacancy defects, the excess (positive) charge on the 
nearest Mo adatoms around the vacancy slightly decreases, and hence does not cause any magnetic 
moment. The charge transfers and resulting excess charges around Mo and M0S2 vacancies are 
affected even up to third nearest neighbor atoms. For example S atoms around Mo vacancy have 
0.1 electrons less charge, since the third Mo atom, which provides excess electrons to S is missing. 
However these S atoms surrounding Mo vacancy receive more charge from nearest two Mo atoms. 
Similar cases also investigated for MoS vacancy defects. Nonetheless, the charge transfers at the 
close proximity of S, S2, Mo and MoS vacancies are not significant as compared with those of the 
perfect structure and consequently do not lead to a magnetic state. Interestingly, the disturbances 
in the charge transfer due to MoS2-triple vacancy are significant and result in the magnetic state at 
the close proximity of the defect. 



In I [figure] [5][]5[ we present the isosurfaces of the difference of charge density of spin-up and 
spin-down states (i.e. Ap^ ^ = p^ — p^ at the close proximity of M0S2 vacancy. After reconstruction 
around the vacancy, 2 Mo and 4 S atoms have dangling bonds. In this case, Mo atoms having dan- 
gling bonds are also less positively charged and S atoms having dangling bonds are less negatively 
charged as compared with those of perfect M0S2. However, in comparison with S, S2, Mo and 
MoS vacancies the amount of charge transfers here are almost doubled to cause to significant dis- 
turbances and spin polarization. The total magnetic moment of 2 /i# are originated equally from d yz 
and orbitals of Mo and p orbitals of S which have dangling bonds as seen in | [figure] [5][]5[ The 
nonmagnetic state is ^ 130 meV energetically less favorable. These results are also consisted with 
the vacancy defects in armchair edged M0S2 nanoribbons. 11 Electronically, vacancy defects give 



rise to states in the band gap, which are localized at atoms around the vacancy (see |[table][3][]3| ). 
The band gap, as well as the electronic properties of IH-M0S2 are modified by these states. 
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Figure 5: Isosurfaces of difference charge density of M0S2 vacancy defect in the (7x7) supercell of 
IH-M0S2. Dashed atoms and bonds are vacant sites. Difference charge density is obtained from 
the difference of spin-up and spin-down charge densities. (Ap^ ^ = p^ — p^) The total magnetic 
moment is calculated as 2 Up arrow indicate the excess spin-up charge. Isosurface value is 
taken as 3x1 0~ 3 eV/A 3 . 

Discussion and Conclusions 

The adsorption of adatoms and creation of vacancy defects in 2D single layer of M0S2 honeycomb 
structure have crucial effects on the electronic and magnetic properties. We found that several 
adatoms can be adsorbed readily at diverse sites with significant binding energy. In this respect, 
M0S2 appears to be a material, which is suitable for functionalization. While two dimensional, sin- 
gle layer M0S2 is a direct band gap semiconductor, magnetic or nonmagnetic localized gap states 
due to adatoms occur in the band gap and expand the capacity of this material in nanoelectronics 
and promise future applications. Additionally, two dimensional sheets of M0S2 can attain a local 
magnetic moments through adsorption of 3d-transition metal atoms, as well as Si and Ge. Car- 
bon being of particular interest, we also examined the adsorption of C2 and linear C3. Significant 
amount of charge is transferred to (or from) adatom. Excess surface charge at higher coverage of 
adatoms can improve tribological and catalytic properties of IH-M0S2. While vacancy defects of 
S, S2, Mo and MoS created in 2D IH-M0S2 do not induce any magnetic moment, the creation of 
M0S2 triple vacancy results in a significant magnetic moment in the system. Briefly, functionaliza- 
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tion of M0S2 honeycomb structure through adatom adsorption and vacancy creation appears to be 
a promising way to extend the applications of M0S2. Since the interlay er interaction in graphitic 
M0S2 is weak van der Waals, present results are also relevant for the surfaces of 3D M0S2. 
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